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Abstract 

The basic principles of the correlation femtoscopy, including its correspondence to the Hanbury 
Brown and Twiss intensity interferometry, are re-examined. The main subject of the paper is 
an analysis of the correlation femtoscopy of small systems with typical sizes about 1 fm. The 

at! 

. uncertainty principle plays an important role for correlation femtoscopy of such systems. It is taken 

into account by the use of the formalism of partially coherent phases in the amplitudes of individual 
emitters. The general analysis is illustrated analytically for the case of the Gaussian approximation 
JL . for emitting sources. The reduction of the interferometry radii and suppression of the Bose-Einstein 

at! 

correlation functions for small sources due to the uncertainty principle are found. There is the 
positive correlation between the source size and the intercept of the correlation function. The 
peculiarities of the non-femtoscopic correlations caused by minijets and fluctuations of the initial 
states of the systems formed in pp and e + e~ collisions are analyzed also. The factorization property 
for the contributions of femtoscopic and non-femtoscopic correlations into complete correlation 
function is found. 
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I. INTRODUCTION 

The correlation femtoscopy, or intensity interferometry method, is the direct tool to 
measure the spatial and temporal scales of extremely small and short-lived systems created 
in particle and nuclear collisions with accuracy of 10 -15 m and 10~ 23 sec correspondingly. 
The method is grounded on the Bose-Einstein (BE) or Fermi-Dirac (FD) symmetric 

properties of the quantum states. It has a deep analogy with the intensity interferometry 
telescope that was proposed by Hanbury Brown and Twiss for measurements of angular 
sizes of remote stars |[4j]. In distinction on standard telescopic and microscopic techniques 
based on the registration of intensities of light or particles, e.g., electrons, coming from (or 
through) the object, this method deals with the correlation between intensities of the source 
radiation registered by two (many) spatially separated parts of devices such as telescopes, 
reflectors, particle detectors, etc. In fact, it measures the correlations between numbers of 
emitted identical particles detected in separated parts of the detector. 

The femtoscopic space-time structure of the systems is typically represented in terms 
of the interferometry radii. They are result of the Gaussian fit of the correlation function 
defined as a ratio of the two- (identical) particle spectra to the product of the single-particle 
ones. In the pioneer papers the measured interferometry radii were interpreted as 



the geometrical sizes of the systems. Later on it was found 



-tZ] that for typical systems 



formed in experiments with heavy ions, the above geometrical interpretation needs to be 
generalized. The treatment of the interferometry radii as the homogeneity lengths s), l| 
in the systems and crucial suggestion as for femtoscopy scanning of the source radiation in 
different momentum bins bring the possibility to analyze different parts of the source and 
explain the behavior of the interferometry radii. In addition, the practical method how to 
use the final state interactions (FSI) and effects of long-lived resonances to extract the B 



correlations in relatively large systems created in heavy ion collisions has been proposed 



io| 



An another challenge, which is still actual, concerns the femtoscopy analysis of relatively 
small systems created in particle interactions such as pp and e + e~ where the observed femto- 
scopic scales are approximately 1 fm or smaller. Typically, the suppression of the correlation 
function is fairly large in these processes. Here we will analyze the femtoscopy of such small 
systems accounting for the uncertainty principle, coherence of the radiation from spatially 
very closely set emitters and non-femtoscopic (non BE, FD and FSI) correlations. The latter 
appear due to the energy-momentum conservation law and incoherent contributions to the 
two- and single-particle spectra induced by particles clusterization in momentum space and 
fluctuations of initial conditions of the collisions processes. The detailed analysis of these 
theoretical problems can help to provide the correct femtoscopy study of the small systems. 



II. THE BASIC IDEAS OF INTENSITY INTERFEROMETRY TELESCOPE. 



The intensity interferometry method for the stars angular sizes measurement was pro- 



posed and realized first by Hanbury Brown and Twiss in [4] at the end of fifties of XX 
century. The electromagnetic radiation from the star is the mixture of different (partially) 
monochromatic wave trains which are mutually incoherent at the moment of radiation. To 
see the principal aspects let us consider the emission from the different sites of a radiating 
object. If the stellar object is close to observer, like our Sun, then one can easily select, 
say, the opposite sites (edges) of it. If the two telescopes are directed to those different 
sites and one measures the correlations between photon numbers coming to each of the two 
telescopes, then the waves from different sites of the source do not mix in the telescopes 
and the correlations between them are absent - the signals coming to the two telescope 
reflectors are mutually incoherent. If, however, the stellar object - let us consider now the 
double star system - is very remote, so that it is impossible to select only one of the two star 
by a telescope, the light from both stars will come to each of the two reflectors and become 
mutually (partially) coherent. For simplicity let us imagine that at some moment of time 
the telescopes register only one wave train from each star, and both these trains are equally 
polarized, and follow each other continuously. Ipso facto we ignore the real problems of the 
intensity interferometry method - how to extract the signal from noise, but preserve the 
principal point of this method. 



A. The basic formalism 



One can decompose the electric field strength into positive and negative frequency parts 
E = EW(x,t) + EH(x,t), E< - )(x,t) = ^^'(x,^. The ideal photon counter reacts just 
to the product E^)(x, t)EW(x, t) = |EW(x,t)| 2 [11]. Far from the stellar object the light 
is described good by the plane waves, so that 

E (+) (x, t) = Ae^ 1 '™*) + Be^^ (1) 

The complex amplitudes A and B have stochastic independent phases: A = \A(t)\e i<PA ^' x \ 
B = \B(t)\e iVB ^ which are roughly constants during coherence time T co h of the wave train, 
e.g. r coh « 10~ 8 sec, so that being averaged over period of time T » r coh the amplitudes 
and their product become zero: (A) = 0, (B) = 0, (AB) = 0. The intensity of light, that 
is proportional to the number of photons registered by the telescope/reflector and photo- 
multipliers at point x, is 



(W*)> = <E 2 > = <E(-)(x,t)E«(x,t)) = <|E(+>(x,t)| 2 ) 

(2) 

= (\A\ 2 + \B\ 2 + ^4_B* e i (( k i _k2 )' x_ ( a;i_u;2 )') -f A* Be'^ 1 ^ 2 ^^^ 1 ^ 2 ^) 
= \A\ 2 + \B\ 2 

where we supposed that \B(t)\ are almost independent on t. 

The statistically averaged intensity registered by one of the telescopes contains the infor- 
mation only about (averaged) squared modulus of the amplitudes. However, the correlation 
of intensities C, defined as the ratio of the averaged product of intensities registered at the 
two space-time points (xi,ti) and (x 2 ,t2), to the product of averaged intensities registered 
at these points, depends already on the differences of momenta and energies of the light 
quanta. To simplify notation we put \A\ = \B\ and get 

c = (/ kl , k2 (xi)/ k k2 (x 2 )) = I ( |^( +)(x t)n ^ +)(xa t) , 3) 
(/ kl ,k 2 (xi)) (/ kl ,k 2 (x 2 )) 4|A| 4 

= 1+2 cos[(ki-k 2 )(xi - x 2 ) - (u! - u 2 )(ti - t 2 )\ (3) 
~ 1 + - cos{6 |k| d + (x L i - x L2 -ti+ t 2 )(wi - uj 2 )\ 



where one took into account that (|A| 2 A*5 = 0), etc. Here the 9 = D/L is the angular 
size of the double star system with "transverse" distance D between stars in the plane 
perpendicular to the direction to the system, L is the distance to the system, d is the 
"transverse" distance between two telescopes, |k| is mean detected wave number, xn is the 
"longitudinal", directed to the system, coordinate of the z-telescope and t\ — t 2 = A is the 
signal delay between points xi and x 2 . It is worthy of noting that the time resolution r res 
in the method has to be smaller than the coherence time, T res < r co h, in order to provide 
correlation measurement of photon numbers during the mutual coherence time r co h of the 
electromagnetic waves in points x x and x 2 . For stationary process the averaging over large 
period of time T plays role of the averaging over the ensemble of events with duration time 
T co h. Also note that the condition of validity of formula ([3]) is 



Otherwise (A(x 1; ti)A{x 2 , ^2)) = 0, (-B(x 1; ti)B(x.2, £2)) = and the correlations disappear, 
C — 1, because of mutual incoherence of waves coming to the two telescopes. 

B. The nature of HBT effect 

The formula (EJ) demonstrates the idea how to measure the differences in momenta (and 
energies) of photons radiated by remote stars by measuring the correlation function de- 
pending on distances between telescopes/reflectors and time delay. In fact, the momentum 
difference in transverse plane is connected with angular size of a stellar object. As for the 
difference in the energy of photons, it is possible, in principle, to measure it with restriction 
given by Eq. (j4]), and this difference would be associated with different temperatures of the 
two stars, if such a situation could take place. However there is no direct connection of the 
difference in energy of radiated photons with the time and space scales of the stellar system 
and, therefore, only angular size of this object can be extracted in this way. The latter is 
the basic application of the intensity interferometry telescope, and the method was used to 
measure the angular sizes of single remote stars 0. 

1 In this case one deals not with two emitters, but with many emitters N at the whole stellar disk, and 
formula ([3]) is modified at t 1 = t 2 as follows: C — 1 oc E^;=i ex P ['(^Ti — kxj)( x Ti ™ XT2)] N —f° 
( 2J gjk|li ' ' Determining the distance corresponding to the first zero of the correlation function in this 
representation Hanbury Brown and Twiss have measured the angular size of Sirius and some other stars. 




(4) 
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Despite the classical description of Hanbury Brown and Twiss (HBT) method, the de- 
tailed analysis of the measurements accounting for the principle of photon registration 
pointed to the quantum nature of the effect The relationship between classical and 

quantum descriptions of the electromagnetic waves has been established through the for- 
malism of coherent states 11] . Appealing to quantum nature of the electromagnetic fields, 
one can say that the method of the measurement of the star angular size is based on the 
positive correlations between numbers of photons registered in two close space-time points 
because of Bose-Einstein statistics for these quanta. 

If a double star system is so far from detectors that the latter can register only a few 
photons per r res , then one has to use the amplitudes of registration at points x\ and 
%2 of the two photons emitted with mutually random phases from the two stars a and 
b: A{xi) = (a\xi) + (b\xi) and A{x 2 ) = (a\x 2 ) + (b\x 2 ). Then taking the ratio of the 
averaged modulus squared of symmetric (over photons permutation) two-photon ampli- 
tude A(xi,X2) = A(xi)A(x2) to the product of the averaged single-photon amplitudes 
one can get again the correlation structure (Ej). Note that describing a registration of 
the two neutrino from the star one should use antisymmetrized amplitude A(xi,X2) = 
(a\xx) (b\x 2 ) - (a\x 2 ) 

III. THE BASIC IDEAS OF THE CORRELATION INTERFEROMETRY 

In particle physics the positive correlations between numbers of identical pions with close 
momenta emitted from an interaction region in proton-antiproton annihilations were found 

n 

in 1960 by Goldhaber et al. (GGLP effect) |l[ . It was understood that the nature of the effect 
lies in quantum statistics for identical particles demanding the symmetrization of bosonic 
wave function. Later on, basing on this fundamental principle, Kopylov and Podgoretsky 
developed the method of pion interferometry microscope, or correlation femtoscopy. They 
found an analogy between the correlation interferometry and Hanbury Brown/Twiss (HBT) 
stellar intensity interferometry Q|. As one can see below, the mathematical structure of 
the principal aspects of the two methods can be presented in identical form. However, 
physically the methods are "orthogonal": HBT method measures differences of momenta 
and energies of photons radiated by stellar objects while the correlation femtoscopy measures 
the space-time separation of the emission points. Correspondingly, the former is based 



on the analysis of particles correlations as they depend on the space-time separation of 
the telescopes/reflectors, while the latter deals with dependence of the correlations on the 
particles momenta differences. In this sense the term "HBT radii" that one often uses to 



present the femtoscopy measurements is not quite correct, as R. Lednicky stresses 



12|. 



A. Standard approach 

The basic ideas of the correlation femtoscopy are described in many publications, e.g., Ref. 



13l | . We reproduce them here with, however, some important remarks. Let us suppose that 
two identical bosons (e.g., pions) are emitted from the two space-time points, X\ = (ti,Xi) 
and x 2 = (£25X2) and then propagate freely. The wave function of a single particle at the 
initial time t; t in the configuration representation is <5 3 (x — Xi). At some time t in the 
momentum p = (p° = E, p)-representation it is ip x .(p,t) = ^ 2 J^ 3/2 e~ lEt e~ t:x - iP . 

In the momentum representation (pi,p2) the two-boson wave function is symmetrized 
and has the form 

^(pi.ftH t) = [e^e*** + e'^V^ 2 ] e*& +B >* (5) 

Then the probability to find the two pions with momenta pi, P2 is expressed through the 
scalar product of differences of 4-momenta and 4-coordinates of the two-pion emission: 



W Xl}X2 (p 1 ,p 2 ) = \^x u x 2 (pi,P2]t)\ 2 oc 1 + cos[(pi -p 2 ) ■ (x x -x 2 )} ■ (6) 

Comparing Eq.flH]) with result (jSJ), related to the interferometry telescope, one can see 
that they differ from each other by the factor 1/2 before the cosine. The attentive reader 
can notice the difference between the two cases. In the case of emission from a double star 
system the correlation of intensities accounts for all the possibilities: this is the correlation 
between intensities of the wave trains coming from different stars as well as from the same 
star. In opposite, when formula ([6]) is derived, it is supposed that if one boson (with 
momentum pi) is emitted from the point xi, then the another boson (with momentum P2) 
is emitted from different point X2 and vice versa, but the possibility when the two bosons 
are emitted from the same point, xi or X2, is excluded. In the case of independent particles 
radiation all the possibilities have to be taken into account. For such a case one usually 
demonstrates the idea of the correlation femtoscopy method by means of factorization of the 
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two-particle normalized emission function, p(xi,x 2 ) = p{xi)p{x 2 ) and integrates the two- 
point probability fl6]) over space-time region. Then for the correlation function C(pi,p2) one 
has (we ignore here possible correlations between coordinates and momenta of the emitted 
particles) 

W(p u p 2 ) 1 



C(puP 2 ) 



W( Pl )W(p 2 ) 
1 

1 + 



W( Pl )W(p 2 

d 4 xp(x)e iqx 



d x x d x 2 p(x 1 )p(x 2 )W XUX2 (pi,p 2 ) 

2 



(7) 



W(p x )W(p 2 

where q = p\ — p 2 . So, the probability to find the particles with momenta pi, p 2 at very 
large times t — > oo is expressed through the Fourier image of the emission function. This is 
the typical basis of the correlation interferometry method allowing one to analyze size and 
shape of small systems. If p is the Gaussian-like emission probability that in some reference 



frame has the form p(x) oc exp 



E 

i=l 



2R? 



5{t — t ), then ([7]) reads as 



C(pi,p 2 ) = 1 + exp 



i=l 



(8) 



For example, the typical resolving width of the correlation function q = 50 MeV corresponds 
to the size 4 • 10 -15 m = 4 Fm (Fermi) = 4 fm (femtometer) . The origination of the name 
of the method - correlation femtoscopy is obvious from such estimates. 



B. Correlation femtoscopy formalism under microscope 

Problem however appears when we apply the basic formula ([7]) to the system with small 
number of emitters. For example, if there are the only two different emitting points: p(x) = 
\(5 4 (x - x x ) + S\x - x 2 )), then from flTJ follows = 1/2) 

"1 



W xltX2 (pi,p 2 ) OC ^ PiPji 1 + cos [(jPi - Pj)( x i - x j)\) = 1 + cos2 



2 (Pl -P2)(xi ~X 2 ) 



(9) 



The result is incorrect as we shall show. 

To analyze the situation in detail let us consider single particle radiation from the two 
points. If the emission from the point x\ is associated with the quantum state that is 
distinguishable from the state corresponding to emission from the point x 2 , then the two 
amplitudes A, are allowed to be realized with some probabilities pi and one can build the 
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corresponding diagonal density matrix: 

p x : A 1 (p)=e ip ^e- iEt and p 2 : A 2 {p) = e~ lEt ■ p 1 + p 2 = l. (10) 

Here and below we omit multiplier (27r) _3//2 since these factors cancel in the correlation 
function (see below). Note that W XljX2 (p) = P iAi{p)A\{p) + p2v4 2 (p)A 2 (p) = 1. In the case 
of the two-boson independent radiation from the two points there are the three different final 
states (amplitudes) A^ that are realized with the corresponding probabilities py-: 

Pu: A ll {p l) p 2 ) = e^e^e'^ +E ^\ p 22 : A 22 ( Pl ,p 2 ) = e ^ e ^ e -^+^)t. 

(11) 

Pi2- A l2 ( Pl ,p 2 ) = 4= (e ipiXl e ip2X2 + e ^ X2 e ip2Xl ) e -^+^)* ; Pn + p 22 + p 12 = 1. 
v2 

If pa = p\ and pi = p 2 = 1/2 then p n = p 2 2 = 1/4 and p% 2 = 1/2. Then the probability 
to find the two particles with momenta pi,p 2 that are emitted independently from the two 
points x% and x 2 is 

C(pi,p 2 ) = u/ Wx )7\w V2 \ n \ = Piil^i| 2 = l + ^ cos [(Pi-P2)(si-a;2)]. (12) 

The last result coincides with Eq. ()3]) for stellar intensity interferometry telescope. The 
result f JT2|) takes into account that the emission of both quanta from the same point (say, xi) 
brings no interference while the wrong account for the interference leads to the formula 0. 
The latter becomes true only for large number of independent emitters. Therefore the 
simplest "derivation" of the intensity interferometry method (J7J) has to be corrected to 
exclude the double accounting in the correlation function for the contribution associated 
with the particles pairs emitted from very close points. This correction is significant if the 
number of independent emitters is not large. We shall consider this effect in detail in the 
next section. 

If the emission from points x\ and x 2 is not independent at all (full coherence), then 
one cannot use the diagonal matrix density associated with the probabilities pi,p2 but only 
the pure state: A xljX2 (p) = Ai(p) + A 2 (p). The amplitude of the two identical bosons 
symmetrized over p u p 2 in this case is: A 3 . UX2 (p 1 ,p 2 ) = A^^p^A^^p^ and 

n( \ W XltX2 {pi,p 2 ) 

C{Pl,P2) = 777 7 Vyi/ FT = L 13 

W XUX2 {pi)W x ^ X2 {p 2 ) 
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C. Correlation femtoscopy in random phase representation 



Both these cases of completely independent and fully coherent radiations from points x\ 
and x 2 can be reproduced in the formalism of partially coherent phases jl^ . To demonstrate 
it let us consider the two-point source with coordinates x± and x 2 and some undetermined 
phases 4>(xi) (i = 1, 2) and express amplitude of single boson emission with momentum p 

A XuX2 (p) = (e tpxi e^ {xi) + e 4 P x V^ 2 ))e-^ (14) 

The symmetrized amplitude of the two-boson radiation is 

A Xl ,x 2 (Pl,P2) = Aci, !B2 (Pl)A B i,x a (P2)- (15) 

The probability to register the two identical particles with momenta pi and P2 is 



W( Pl ,p 2 ) =< A^MA^MAl^ip^Al^pz) > (16) 

where brackets mean the averaging over events with partially coherent phases (over all events 
of the emission). Then in the case of random phases, ^e^* 2 '^ 1 ^) = 5 4 (xi — x 2 ), we get 
result f ll2p and for fully coherent radiation when /e l ^ X2 '~^ Xl 'n = 1 we get f )13p . 

In the case of independent emitters in points X\,x 2 the above phase average can be 



presented in the form 14j: 



G 12 = ( e ^M-^i))) = 8\x x - x 2 ) = I 12 5(h - t 2 ) (17) 
where l\ 2 is the overlap integral 



1*3 



rf 3 x^(t,x)V*.(t,x) 



(18) 



As it is demonstrated, an independent emission with fully random phases corresponds to the 
mixed state with probabilities pi for z-states. Such a description is possible only if bosons 
are emitted independently from different points Xi in distinguishable/ orthogonal quantum 
states. The latter requirement is satisfied in above case of the flat momentum spectrum for 
each emitter, f(p) = W(p) = const, since the initial quantum states <5 3 (x — xi) taken in 
different points i are orthogonal. 
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D. Uncertainty principle and formalism of partially coherent phases 



If the momentum spectrum f(p) is essentially not flat, this corresponds to the wave packet 
characterized by its center X{ and some finite width. One can discriminate between the 
different states % and j only if they are approximately orthogonal: the overlap integral ( TT8|) 
is small, I{j <C 1. In other words, the distance between the centers of emitters has to be larger 
than the width of the emitted wave packets. Since the latter is the inverse of the variance 
Ap of the momentum spectra, so (xj — Xj) 2 > 1/Ap 2 . The latter expresses the uncertainty 
principle: one can discriminate the wave packet % without noticeable violation of the particle 
spectrum, if the measurement, that localizes the particle's position somewhere inside the 
sphere with the center x^ and the diameter not less than 1/Ap, unambiguously points to 
the quantum state i, but not j. So, the distance Ax = (xj — Xj) should satisfy uncertainty 
principle in the form Ax 2 Ap 2 > 1. Then the system presented by the superposition of 
the quantum states Ai with phases </>(xj) like in Eq. (|T4l) . after a phase averaging ( IT71) is 
approximated well by the mixed state with diagonal density matrix py. 

In relativistic physics there is another uncertainty principle: the measurement of the 
particle's momentum p has accuracy depending on the duration of the measurement 5t: 
5p ~ l/5t [l6|. So one can measure the time of particle emission without noticeable violation 
of the momentum spectrum with accuracy not better than 1/Ap and, so, the ^-function in 
Eq. ffTTj) has to be smeared when one deals with the wave packets. So, for normalized wave 
packets, one cannot use the random phase approximation if the distance in space and time 
between emitters is less than the width 1/Ap of the wave packet (in units h—c—1). For 

1 - p2 Ap 2 (x 1 -x 2 ) 2 

example, if /(p) = (27rAp2)3/2 e W , then J 12 = e 2 at tx = t 2 and 5(ti - t 2 ) =>• 

_Ap 2 (i 1 -i 2 ) 2 

^12 — ^ e 5 • The last term expresses uncertainty principle for momentum-time 
measurements. 

Therefore, the fully chaotic phases, or diagonal density matrix, are possible only under 
some conditions which we discussed before. For normalized quantum states emitted from 
the points i and j the phase average can be written as 

Gij = (e^**)-^))) = I ije - AP (t i~ tj) (19) 

where is the overlap integral (11 8p . and Ap is a variance of the momentum spectrum of 
emitters. 
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In opposite case when (xi — Xj) 2 <C 1/Ap 2 the states are indistinguishable and overlap 
integral f lT9|) ~ 1 at t\ ~ t?, since at xi = Xj it is just normalization. Then we come back 
to the fully coherent radiation ( I13p which takes place for very close emitters 15j. So, in 
both limited cases of chaotic and fully coherent emission Eq. ( JT§j) leads to physically obvious 
results, and so we will use Eq. (|T9|) at any distances between the emitters as accounting for 
uncertainty principle. In this way we solve the problem of distinguishability (or indistin- 
guishability) of the quantum states ipi associated with different emission points xf. we always 
sum i-amplitudes, not i-probabilities, and average the resulting distributions/spectra over 
partially coherent phases in correspondence with the requirement of maximal possible dis- 
tinguishability and independence of different z-states compatible with uncertainty principle 
for momentum & position and energy-momentum & time measurements. 



E. Femtoscopic correlations in multiparticle systems 

We discussed above the simplest situation with the only two emitted pions. Typically, the 
basic formalism of correlation interferometry in multipion systems for inclusive momentum 



spectra is similar 17] : one can use the same basic formula (j7J) with the substitutions: 
W(p) -> p°f£, W( Pl ,p 2 ) -> rf^^fe, p{x) -> S(x,p = ^). The correlation between 
the space-time position x of the emission point and momentum p of emitted particle, reflected 
in p-dependence of the emission function S(x,p), leads to dependence of the interferometry 
radii in formula (jSJ) on the mean momentum p of pion pair, Ri — > Ri(p) |6|]. The x — p 
correlation appears due to fast expansion of the multiparticle systems created in high energy 
A+A (or even p+p) collisions. The radii Ri(p) are associated then with homogeneity lengths 
in i-direction in inhomogeneous expanding systems [8|. So, one can carry the two-pion 
correlation results from the two-pion system to multiparticle systems, if applies the above 
substitutions. In fact, our quantum- mechanical consideration is related to the rest frame of 
the homogeneous emitting subsystem forming the spectra in a vicinity of some momentum 
P- 

The another aspect of multiparticle emission is that the two-pion wave function may 
not be always factorized as a subsystem: ^ = |7r + 7r + ) \X), where state X is related to the 
residual part of the total system. The symmetrization/antisymmetrization procedure has to 
be applied to the total system wave function However, it can be difficult to provide it 
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if one supposes the distinguishability of the radiation points as in the standard correlation 



interferometry method of independent sources l[ 



2|. For example, let us suppose that one 



pion with momentum pi, Tt + (p±), is the primary particle emitted from the fireball near the 
point X\ and another pion 7r + (p 2 ) with the momentum p 2 is emitted together with particle 
X at the resonance decay, Res — > vr + (p 2 ) + X(p 3 ) (say, p° — > 7r + (p 2 ) + n~(P3)) near the 
point X2- Such a system may not have 7r(pi) <-> 7r(p 2 ) symmetry because 

i) if at the vr(pi) <H- 7r(p 2 ) exchange the momentum p 3 of the particle X is preserved, 
then the momentum of emitted resonance is changed and gets the other value pi+p^. This 
new resonance state has typically another probability to form as compared to the previous 
resonance state. (In extreme case it can even happen that the value (pi + p%) 2 instead of 
v\es = G°2 +P3) 2 , excludes the resonance decay into the pair ir + { m ) , X(p 3 ) because of the 
kinematics). In general all that suppresses the interference effect [181 ]; 

ii) if the exchange of tc + momenta is accompanied by the change of p 3 by the value 
q = Pi — P2, Ps — > P3 = P3 + q, then the two configurations |vr+(p 1 )7r+ (p 2 )X(p 3 )) and 
\TT^ l (p 2 )7r^ 2 (pi)X (p' 3 )) become distinguishable and so cannot interfere; 

iii) if at tt(pi) o 7i(p2) the particles X(p 3 ) together with ir + (p2) are emitted by the res- 
onance Res from the fireball near point X\ and primary pion vr + (pi) radiates from x 2 , the 
requirement of distinguishability of the emission points may lead to principal distinguisha- 
bility of the points of emission X\ and x 2 of the particle X that excludes the pion interference 
effect. 

So, the exact symmetry of the total system at p\ -H- p 2 exchange can be lost in similar 
situations and the two configurations ^(pi)^ (p 2 )) and [^(/^^(Pi)) of the two-pion 
subsystem can interfere only partially. This has to be taken into account in the transport 
models used for high energy collisions. 

The attempt to estimate the limits of applicability of the standard correlation femtoscopy 
is done in Ref. Ij]]. It is argued that when W^iP^^iP^)) Wxi (P2)tt+ (pi)) then the 
classical phase-space position of the residual part of the system is changed by the value 
q(x\ — x 2 ) along each direction i. If this change exceeds the size 2irh of the elementary cell 
of the phase-space per one degree of freedom, then the residual system moves to another 
quantum state, and it can be, in principle, measured. Therefore the pion interference should 
disappear when qiRi > 2tt (H = c = 1), where Ri is the effective system size /homogeneity 
length in z-direction. However this argumentation fails already in the simplest case when 
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the system contains only identical particles. Such a system can be symmetrized with distin- 
guishable and independent radiation points, and there is no principal possibility to measure 
the exchange of the momenta between two identical pions at any value of q using for this aim 
the residual system containing the same sort of pions. We think that the real picture of the 
pion interference can be restored on the basis of symmetrized/antisymmetrized amplitude 
of the total system with partially coherent phases in the way described in subsection IID. 



IV. THE CORRELATION FEMTOSCOPY IN GAUSSIAN APPROXIMATION 

Here we apply the basic ideas discussed in the previous section to construct the simple 
analytical model accounting for uncertainty principle in the correlation femtoscopy. We 
use the non-relativistic approximation in the rest frame of the source which has Gaussian 
sizes corresponding to the homogeneity lengths Ri in the corresponding part of the total 
expanding system. The Lorentz transformation of the results to a global system where 
source moves with 4-velocity is trivial enough. Let us introduce a single bosonic state 
ip x (p, t) corresponding to emission of a single particle with mass m and 4-momentum p from 
4-point x and its free propagation: 

^CM) = e* (px - p0 V^/(p) (20) 

where tp(x) is some phase and / defines the primary momentum spectrum /(p) that we take 
in the Gaussian form with variance Ap = k: 

/(P) = / 2 (P) = (d^S«-*- (2D 

The amplitude of the single-particle radiation from some 4-volume at very large times 
t^ can be written as a superposition of the wave functions ip x {p) with some coeffi- 
cients/distribution p(x) for emission 4-points x. In the momentum representation it looks 
like 

A(p,t) = c J d 4 x^ x (p,t)p(x) (22) 

where c is the normalization constant. Let us select the two directions: parallel to beam 
axis, it is marked as index "L" and orthogonal to it - transverse axis "T". The distribution 
function of the emission centers is supposed to be the Gaussian one with normalized to unity 
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quantity p 2 (x, t) which is the probability distribution in the case of random phases (j>(x). 

" w = ^kr (23) 

The single-particle momentum density averaged over events with different phase distribu- 
tions is 

W\p) = c 2 J d A xd 4 x'e ip{x - x ' ) p(x)p(x')(e iMx) -^ x ' ) )f(p). (24) 

To calculate the phase average one needs first to calculate the overlap integral f fl8|) . In non- 
relativistic approach the wave function of the particle emitted from the point Xj at moment 
ti in coordinate representation is at some moment t 

^(t,r) = j^j- 2 J f(p)e-^-^e-^^d 3 p. (25) 

Then the modulus of the overlap integral ffl8|) is 



fc^x.-x^.) 2 



dh^ Xt (t,r)r x At,r) 



g 2(l + fc4(t i -t J -) 2 /m^) 



(26) 



To provide calculations in analytical form we substitute the squared time difference 
(ti — tj) 2 by the constant proportional to its mean value over emission region 

(U - tj) 2 a((ti - tj) 2 ) = a J d A x i d A x j p(x i )p(x j )(t i - tj) 2 = AaT 2 = aT 2 (27) 

As we will demonstrate later, the results obtained within such a prescription reproduce 
with good accuracy the exact numerical calculations for momentum spectra and correlation 
functions with the values of a depending on the basic parameters of emission. Then the 
overlap integral f fl9|) that accounts for uncertainty principle takes the form 

- fc 2 (x-x') 2 

G xx , = <e^)-^'')> = pgg^g „-tf(t-n'/2 (2g) 

yl + ak 4 T 2 /m 2 

Now one can obtain the one-particle spectrum (|24|) . It is presented below for the case 
when homogeneity lengths are equal, Rt = Rl = R' 

_ p 2 2p 2 fl 2 p 4 T 2 

W(p) = Ne i+^o* 2 2m2 ( 1+4fe2T2 ), (29) 

where 

k 2 = k 2 /(l + ak A T 2 /m 2 ). (30) 
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The two-particle spectrum averaged over events with partially coherent phases is 
W( Pl ,p 2 ) = c 4 J d^AdV^V^ 1 * 1 ^ ■ 

The 4-points phase correlator supposing emitters to be chaotic and independent in a maxi- 
mum possible way permitted by uncertainty principle is decomposed into the sum of products 
of the two-point correlators (JTHJ) and contains the three terms 

(ei (,(,0+.(, 2 )-,W)-.(4))) = (;,.,,..(;,.,,.. + G X1X ,G X2 , A - G X1X ,G X2X ,G X1X2 (32) 

where the third term removes the double accounting which appears in the sum of the first 
two terms when x\ ~ x 2 ~ x[ ~ x' 2 . See discussion in the Subsection IIB. 

Let us split transverse direction " T" into the two ones: out that is directed along the total 
transverse momentum of the pair pit + p 2 r and side that is orthogonal to the o^-direction 
as well as to the longitudinal "L" axis. Then the calculations of the one- and two-particle 
spectra lead to the following correlation function in the variables of the half of the bosonic 
pair momenta sum p = Si±£l anc l the momenta difference q = Pi — P2 



W On r> ) „2 ff 2 Ak l R T „2 R 2 Ak l R L (qp) 2 T 2 4fc 2 T 2 

' ' IP, q) — — 1 + e i+4fc « T 1 + 4fc "i m — C d (p, q), (33) 

^(^1)^(^2) 

where W(pi) is defined by (p4l) . fc is ([30]) and subtracted function eliminates the double 
accounting at the averaging of the 4-points phase correlator, it is associated with the third 
term in Eq. ( 1321) . This term reduces the intercept of the correlation function as one can see 
below. 

The coefficient a in the expression for fc should be chosen from the requirement of a 
good agreement between the correlation function calculated using the approximation ( 128|) 
and the one calculated numerically using the exact expression for the space part of the phase 
correlator ( )26l) . For the typical freeze-out temperature Ty . = m n = m the value of k is 
k = ^Tf M jn = 0.14 GeV. Having fixed the value of this parameter we present in Figfj] 
the comparison of the sz<ie-projections of the correlation functions calculated numerically 
using ( 126|) without subtraction of the double accounting and corresponding analytical ex- 
pressions (the first two terms in ( 133]) ) . The values of a at different system sizes are presented. 

16 




FIG. 1. The comparison of the correlation functions in side direction without subtraction of the 
double accounting with corresponding analytical approximations. The a parameter values, which 
give a good agreement with the exact results, are presented for different system sizes. For R=1.5 
and 5 fm the curves are merged, while for R=0.5 fm there is a disagreement with exact result at 
large q. The momentum dispersion k = m = 0.14 GeV, p=0, T=R. 

For fairly small systems, R ~ 0.5 — 1.5 fm they are about unity, a ~ 1 and are decreasing, 
a — > 0, when the system size R grows up. 

As it follows from Eq. fl33|) the observed Gaussian interferometry radii of the system 
are reduced as compared to the standard results (we mark it by st index) for the interfer- 
ometry radii for the Gaussian source. The reductions for side- (S), out- (O) and long- (L) 
interferometry radii in the LCMS system (pl = 0) are 

R 2 S 4k 2 R 2 T 
R% t ~ 1 + 4k 2 R 2 T 

Rl ±klR\ 
Rht i + 4fc 2 i?i 

where v out = p out /m < 1. 

By means of Eq. ( 131~]) one can calculate the subtracting correlation function Cd corre- 
sponding to the third term in Eq. fl32|) . It brings the result 
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2q 2 fc 2 | fl 4 (l+8fc 2 l .R 2 J 



C d {p, q) = F(k 2 R 2 } k 2 T 2 )e ^^W^TE^) ^ 



F(k 2 R 2 k*T*) = /° (l + 4^)^(l + 4^f 2 

k (l + 8k 2 T 2 + 8fc 4 T 4 )V2 + + gk^f 2 



One can note that the pre-exponential coefficient F(k 2 R 2 , k 2 T 2 ) tends to zero at large 
kyR 2 and/or k 2 T 2 and tends to unity when both of these quantities tend to zero. In the 
former case of large sizes and/or hard radiation the intercept of the correlation function tends 
to 2 and in the latter case, when the sizes of the source are very small and/or radiation is 
soft, the intercept tends to 1. For example, at T=0 



k 2 R 2 > 1, C(p, q) = 1 + e~ q R (37) 

which is associated with the standard results, and 

k 2 R 2 <^l, C(p,q))«l (38) 

The last result corresponds to undistinguishable positions of the emitters due to uncertainty 
principle in case of very small sources and/or soft particle radiation. 

We illustrate realistic cases in Fig. [2] for sirfe-projection of the correlation function. De- 
spite the correlation function at small k 2 R 2 looses the Gaussian form presented by the second 
term in (133]) . the Gaussian fits still follow the tendencies of the analytic results (13"4"]) which 
demonstrate the reduction of the interferometry radii when the sizes of the sources become 
small. The elimination of the double accounting leads also to the reduction of the intercept 
of the correlation function when the system size shrinks. So, there is the positive correlation 
between observed interferometry radii and the intercept when the system size is changing. 
The analytic results reproduce the exact ones at a — 2 (R — 0.1 fm), a = 1.6 (R = 0.5 fm), 
a = 0.8 (R = 1.5 fm). It is worthy of noting that at large enough system sizes, R > 3 — 4 
fm, the effects practically disappear. 

V. NON-FEMTOSCOPIC CORRELATIONS 

Considering the correlation interferometry of small systems we cannot pass by the phe- 
nomenon of the so-called non-femtoscopic correlations. For a relatively small systems the 
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FIG. 2. The behavior of the two-particle Bose-Einstein correlation function (side-projection) where 
the uncertainty principle and correction for double accounting are utilized. The momentum dis- 
persion k = m = 0.14 GeV, p=0, T=R. 



additional two-particle correlations may appear in the kinematic region typical for quantum 
statistical (QS) and FSI correlations. The well known example of such additional correla- 
tions is the correlation induced by the total energy and momentum conservation laws (see, 
e.g., Ref. 20|). As opposed to the QS and FSI correlations, which are familiar from the 
correlation femtoscopy method, and, so, are called sometimes as femtoscopy correlations, 
these correlations are not directly related to the spatio-temporal scales of the emitter, and 
are called, therefore, non-femtoscopic correlations. Since the latter affect noticeably the 
correlation functions for small systems, the interferometry radii extracted from complete 
correlation function in e. g. p+p collisions depend on the assumption about the so-called 
correlation baseline - the strength and momentum dependence of the non-femtoscopic cor- 
relations 2ll~l23|. It has an influence on an interpretation of the momentum dependence of 



the interferometry radii in p+p collisions, where the possibility of hydrodynamic behavior 
of matter is questionable. Therefore for successful and unambiguous applications of the 
correlation femtoscopy method to elementary particle collisions one needs to know mecha- 
nisms of the non-femtoscopic correlations to separate the femtoscopic and non-femtoscopic 
correlations. The very important question is whether it is possible to factorize both these 
contributions to the the total correlation function. 
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A. Non-femtoscopic correlations from minijets 



To study this problem we follow the ideas of the Ref. [24| and consider simple analyti- 
cal model of three-particle emission accounting for non-femtoscopic correlations caused by 
minijets fragmentation and energy-momentum conservation law. Even if one ignores the 
QS correlations, the amplitude of the three-particle emission, accounting for minijet frag- 
mentation, cannot be expressed by the simple product of single-particle amplitudes, but is 
modified and can be written as follows 24] 



MPi, P 2 , p 3 ) = A( Pl )A(p 2 )A(p 3 ) -> A( Pl )A(p 2 ) A(p 3 )Q(p 2 , p 3 ) (39) 

where factor Q(Pj,Pj) = exp ^— ^ P ' 7 ^ \ describes the correlation between the particles 
with momenta Pi and p ■ induced by the fact they came from the same minijet (or cluster 
in momentum space). Taking into account different possible cluster structures that can be 
formed by different pairs of particles one can write the single- and two-particle emission 
probabilities in the following way 

W( Pl ) = J2l J dpW2dpl5( Pl -p*)\A(p* l7 p* 2 ,pl)Q(p* 2 ,p* 3 )\ 2 (40) 

w( Pl ,p 2 ) = j2^ J ^^^(pi-pn^-P^i^P^P^pDQCp^pDi 2 (41) 

i¥=3 

And the corresponding correlation function is 

C(p, q) = \C™*{ P , q) + ~C^(p, q) (42) 

where C™ a: (p,q) denotes the contribution to the correlation function corresponding to the 
pairs formed by the "independent" particle and the particle from the cluster while C^ et { P , q) 
is another contribution related to the pairs where both particles belong to the cluster. 

To describe the correlations induced by the momentum conservation law one should 
include certain delta-functions to the expressions for W(p) and W(p 1 , p 2 ). Here we use the 
"soft" form of the momentum conservation law where the delta- function is substituted by 
its gaussian-like approximation, keeping in mind that the conservation law being strictly 
kept for the whole system of produced particles can be only approximately fulfilled for the 
subsystem of identical bosons under consideration. 

S(pl + pj + P3) A(pJ, pS) = Ce-^^ 2 ^ (43) 
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FIG. 3. The correlation function C(q) accounting for the non-femtoscopic correlations due to 
minijet /cluster formation and momentum conservation law, C c i(q), is compared with the pure QS 
correlation function Cbe(q)- One can see that the full CF is quite well factorized into the QS and 
non-femtoscopic parts. The model parameters are 7 = 0.5 GeV, p = 0.35 GeV/c, k = 0.1 GeV/c, 
R = 1 fm 

Fig. [3] illustrates the relation between the femtoscopic and non-femtoscopic correlations 
in our simple model. One can see that for small systems with the size ~ 1 fm the full 
correlation function can be factorized in the two parts corresponding to the femtoscopic and 
non-femtoscopic correlations. 

B. Non-femtoscopic correlations from fluctuating initial state 

The another type of correlations not induced by the quantum statistics effects are the cor- 
relations connected with existing of subensembles of events with different emission functions 
that leads to the corresponding fluctuations in single-particle and two-particle momentum 
spectra. In hydrodynamical models of nucleus-nucleus and proton-proton collisions these 
fluctuations can be caused by asymmetrically fluctuating initial densities used for the hydro 
stage of the model. 

Let us consider the effect of such correlations on the resulting correlation function again 
on the example of simple analytical model. The probabilities of the single- and two-particle 
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FIG. 4. The non-femtoscopic correlation function ("baseline") at the fluctuating initial conditions 
for hydro dynamic expansion of the system. Depending on the relation between value of the pa- 
rameter a, characterizing the strength of the fluctuations and the parameter d = 1 GeV, that 
describes the softening of the momentum conservation law for the pion subsystem, the baselines 
can be growing, constant or decreasing with momentum difference q. 
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FIG. 5. Factorization of the femtoscopic and non-femtoscopic correlations when the latter are 
induced by fluctuating initial conditions for hydrodynamics 
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emission in case of fluctuating emission function should be averaged over the ensemble of 
states corresponding to the events with different initial conditions 

W(jp) = Y i p(u i )W(jp;u i ) (44) 

i 

W(p 1 ,p 2 ) = ^2p(ui)Wi(pi,p2]Ui), (45) 

i 

where Ui denotes the z^-type of initial conditions and p(ui) is the distribution over initial 
conditions, YliP( u i) = 1- ^he corresponding expression for the non-femtoscopic correlation 
function takes the form 

r u „ \ - E i p( u iW(pi;u l )W(p 2 ;u l ) 

Cnf{pi >™ ~ EM^w^EjpMw^uj) (46) 

Let us consider the simple model where the single-particle Mj-dependent distribution 

_ p 2 (p-u) 2 

W(pi,Ui) are obtained using the modified momentum spectra /(p) k e a 1 -> e ^ , 



where p{Uj) distribution is supposed to be the Gaussian p(u) = -^j2 e u2a2 . The momentum 
conservation law we use, as in the previous subsection, in the soft form using the factor 
A(pi,pa) = Ce-( pi+ P 2+ -) 2 / d2 . 

Depending on the ratio between parameters d describing momentum conservation and a 
describing initial conditions fluctuations the baseline formed by the C^p function can be 
growing, constant or decreasing with q (see Fig. H]). 

We see that fluctuating initial conditions for hydro can cause non-femtoscopic correlations 
similar to those resulting from minijet fragmentation. Fig. [5] illustrates the factorization of 
its contribution to the full correlation function. 



VI. CONCLUSIONS 



We discussed the principal problems of the correlation interferometry of small sources, 
in particular, the similarity and the principal difference between Hanbury Brown & Twiss 
intensity interferometry and pion correlation interferometry methods. The main subject of 
the paper is the effects caused by not complete distinguishability of the emitting points in 
the source because of the uncertainty principle. To account for the latter the formalism 
of partially coherent phases in the amplitudes of individual emitters is developed. The 
specific treatment is required for elimination of the double accounting to the correlation 
function from the configurations in four-points phase correlator when all individual emission 
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centers are close. It is shown that the effects are significant for the systems with the sizes 
and emission duration times about 1 fm, and they are expressed in the reduction of the 
interferometry radii and suppression of the Bose-Einstein correlation functions. There is the 
positive correlation between the source size and the intercept of the correlation function. At 
the typical slopes of the momentum spectra in p + p and A + A collisions these effects are 
negligible for fairly large sources with radii about or more than 3-4 fm. 

The peculiarities of the non-femtoscopic correlations caused by minijets and fluctuations 
of the initial states of the systems formed in pp and e + e~ collisions are analyzed also. The 
factorization property for the contributions of femtoscopic and non-femtoscopic correlations 
into complete correlation function is demonstrated in the simple analytic models. 

The application of the theoretical results to pp and e + e~ collisions will be done in the 
separate work. 
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